Herzfeld versus Mott transition in metal-ammonia solutions. 
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Although most metal-insulator transitions in doped insulators are generally vie- 
wed as Mott transitions, some systems seem to deviate from this scenario. Alkali 
metal-ammonia solutions are a brilliant example of that. They reveal a phase sepa- 
ration in the range of metal concentrations where a metal-insulator transition occurs. 
Using a mean spherical approximation for quantum polarizable fluids, we argue that 
the origin of the metal-insulator transition in such a system is likely similar to that 
proposed by Herzfeld a long time ago, namely, due to fluctuations of solvated elec- 
trons. We also show how the phase separation may appear : the Herzfeld instability 
of the insulator occurs at a concentration for which the metallic phase is also uns- 
table. As a consequence, the Mott transition cannot occur at low temperatures. The 
proposed scenario may provide a new insight into the metal-insulator transition in 
condensed-matter physics. 
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In the first theory of metallization, Herzfeld jl| has considered dipolar fluctuations of 
neutral atoms as the origin of a metal-insulator transition (MIT). He emphasized that the 
restoring force of an electron bound to an atom collapses and the substance becomes metallic 
at increased densities due to local field effects. Could the Herzfeld idea be applied to describe 
the MIT in real systems ? Despite some successful examples of such applications 
serious doubts arose because most of the experiments gave evidence in favor of the Mott 
scenario [6j. The latter focused on the screening of the long-range Coulomb potential by a 
stable electron gas, which prevents the formation of bound electrons. As far as this metal 
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Fig. 1: Phase diagram of Na-NFhj. The experimental data on the locus of the phase separation 
are indicated by the square [9( and diamond symbols lOJ, respectively. The triangles show the 
change in sign of the derivative of the conductivity coefficient da/dT, which is used to estimate 



the locus of the MIT Jj|. The solid curve corresponds to our calculations of the locus of the 
polarization catastrophe. All dotted curves are guides for the eyes, except the dashed horizontal 
line at T = —80° C indicating the solidification temperature of ammonia. [l(MPM) ~ 2 • 10 20 
cm -3 ]. 

hallmark process runs, the system remains metallic, but when it ceases at lower densities, 
a MIT occurs. It is very important to realize that the two theories are based on different 
effects, since Mott considered the MIT from the metallic side of the transition, whereas 
Herzfeld investigated the same phenomenon from the insulating one. No comprehensive 
connection between them has been provided up to now. Recent studies on a Wigner crystal 



formed by 
was shown 



arge polarons have however opened new perspectives for the Herzfeld idea. It 
3,ls| that this Wigner crystal of polarons loses its stability owing to the dipolar 
interactions between polarons, yielding a polarization catastrophe, which provokes the onset 
of metallization. This phenomenon is a quantum version of the classical Herzfeld scenario. 
Motivated by these results, we argue that a modified Herzfeld approach may provide a key 
to understand the MIT in certain real systems. In particular, we focus in this paper on the 
metal-ammonia solutions (MAS). 

n 

Although it is a century-old problem [9(, the phase diagram of MAS (see Fig. 1) has 
remained mysterious up to now. Many studies have been performed and a large volume 
of experimental data have been accumulated about this fascinating system (for review, see 
[ill, Q])- Once an alkali metal is dissolved in liquid ammonia, it immediately dissociates to 
give two separated entities with unlike charges : the solvated ions and the excess electrons. 



At low metal concentration, the solution remains non metallic (electro 



ytic) and has an 



intense blue colour independently on the type of alkali metal. Jortner [l4] argued that due 
to short-range interactions with ammonia molecules, an excess electron forms a cavity free 
of solvent in which it localizes with the help of the polarization carried by the surrounding 
ammonia molecules. This process results in a trap formation similar to that for the polarons 

i n 

in solids. The radius of the cavity has been estimated to be r c « 3.2 A [14J. Modern theories 
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based on path integral simulations 15], or on the density functional approach [16|,ll7| provide 
an evaluation of the microscopic structure around solvated electrons but they all yield the 
same physical picture as that described above. At large enough metal concentration, the 
MAS becomes a liquid metal with a typical bronze coloration. However, at concentrations 
varying from 1 to 10 mole percent of metal (MPM), a separation between the low density 
blue phase and the higher density bronze one takes place, resulting in a miscibility gap 
below a critical temperature (Fig. 1). Importantly, the phase separation occurs for Li, Na, 
or K, but was not observed in the case of Cs. However, for all type of alkali metal, many 
experimental data reveal the presence of a MIT in the same range of densities This is 
reported in Fig. 1 for the case of Na. 

n 

Earlier models considering the Mott mechanism [13J] or involving an association of loca- 



lized electrons in clusters 



18] were not able to explain the whole phase diagram observed 



in MAS satisfactorily. What is the reason? Let's give an outlook on the complexity of the 
problem. From an electrostatic point of view, solvated electrons behave more or less like 
some solvated anions, the counterpart of the solvated metal cations. The Debye screening 
length is found to be about A at 4 MPM (taking into account the static dielectric constant 
e s ~ 20 of ammonia), which makes the MAS a strong electrolyte in this concentration range. 
The static Coulomb interactions are thus essentially already screened when the MIT occurs, 
and therefore cannot be its origin. The short-range interactions between electrons are also 
unlikely to be responsible for the MIT and the phase separation, because the mean dis- 
tance between electrons is still about 12 A at the relevant concentration 4 MPM, which is 
enough to neglect any overlapping between the wave-functions of the electrons localized in 
their ground state. Finally, the occurrence of a phase separation at low temperatures in the 
concentration range where the MIT occurs, gives serious doubts about the possibility of a 
Mott transition. We will come back to this point. 

The origin of the MIT must be found elsewhere, and a reasonable hypothesis is that 
it results from quantum momentum fluctuations of the solvated electrons. These ones are 
self-trapped quantum particles, whose dipolar momentum effectively fluctuates due to their 
quantum nature, with a characteristic frequency cjq(T). This frequency corresponds to elec- 
tronic transitions of the electrons between two states bound in their own trap potentials. 
They are experimentally detected by ordinary optical absorption measurements. The lat- 
ter reveal a broad absorption line peaked at ujq{T) ~ 0.9 eV at low concentration and 
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Fig. 2: The density of state (DOS) of the collective polarization modes of interacting dipoles. The 
DOS is drawn in arbitrary units, for various metal concentrations at T = — 35°C. The lower edge 
UJ-(T, n) drives the stability of the system (see text). 

Fig. 3: Concentration dependencies of the dielectric constant (a) and the locus of the maximum 
of optical absorption (b) in MAS. In a) the square symbols correspond to the experimental data 
[27I on dielectric constant at T = +20°C, the circle symbols to that at T = — 35°C 28]. The solid 



and dashed curves show our results at T = +20°C and — 35°C, respectively. In b) the triangle 
symbols indicate the experimental data on the absorption maximum in Na-NH3 at T = — 65°C 
obtained from 29| and square symbols from [30|], while the solid curve shows our results at the 



same temperature. 

T = — 70°C, with a tail extending in the visible region, and providing the blue colour of 
diluted MAS. For a given temperature, the maximum of optical absorption shows a pronoun- 
ced red shift at metal concentration above 0.1 MPM, which cannot be caused by short-range 
or static Coulomb interactions as it has been discussed above. The frequency Uq(T) charac- 
terizing the solvated electron state, is a significant phenomelogical parameter of our theory. 
It is associated to the static polarizability «o(0) = e 2 /muo 2 ~ 10 -22 cm 3 of a solvated electron 
(m is the electron mass). That is a huge pola rizability with respect to the one of a single 
ammonia molecule cxnh 3 ~ 2.8 • 10~ 24 cm 3 [19i] or that of a sodium ion a^ a + ~ 2 • 10 -25 



cm 
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20| | . Therefore, these quantum fluctuations of isolated solvated electrons and their induced 
dipole-dipole interactions may have a dominant role in the MIT. Moreover, since u; (T) is 
much higher than the inverse relaxation time r _1 of the solvent (typically in the THz range), 
the dipolar interactions between solvated electrons are only screened by the high frequency 
dielectric constant of ammonia. 

The key idea of Herzfeld was to evaluate the effect of the local field, with the help of the 
Clausius-Mossoti relation, as a function of the density of a substance and the polarizability 
of its constituents. Nevertheless, his calculation did not take into account the fact that the 
particles interact. However, these interactions induce collective modes of polarization, which 
substantially modify the generalized susceptibility x(u;) with respect to the non-interacting 
case. This susceptibility gives the response to an external field E , i.e. P = x(u;)E , where 
P is the polarization. Taking into account both the local field and the interactions effects, 
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we generalize the Clausius-Mossoti relation as 

e(T,id)/e NHs (T,u)-l 4vr 



XM, (1) 



e(T,w)/e NHa (T,w) + 2 3 
where eNH 3 (T,id) is the temperature- and frequency- dependent dielectric function of pure 
ammonia, e(T, u>) is the similar quantity of the solution. The susceptibility may be expressed 
in terms of an effective dynamical polarizability a (a;) of a single solvated electron, by the 
relation = na(id) / e^- Hence, the problem focuses on the calculations of a(u). The 

speetra, density-funetional theory fl seems to provide a rigorous basis tor sueh ealeulations. 
However, despite recent progresses in this direction [22J, the complexity of the microscopic 
structure and the presence of disorder in MAS prevent accurate numerical calculations of 
the electronic properties. 

To overcome this barrier we choose a simple semi-analytical model suitable to calcu- 
late a(u) with a reasonable accuracy. We consider MAS as a fluid of quantum particles 
localized in cavities with diameter a = 2r c , which interact through induced dipole-dipole 
interactions. They can be treated as a set of quantum Drude oscillators with isolated pola- 
rizability ao(u) — e 2 /m(td 2 — uo 2 ). In the simplest approximation the interactions between 
particles can be cut off at small distances by the cavity size, whereas short-range details 
are ignored. Similar models of quantum polarizable fluids have been extensively studied 
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261 ]. and we use the results which were previously obtained. In particular (see 



Methods), the problem is reduced to solve a quadratic equation, whose roots are complex, 
i.e. a (u) = a'(td) +ia"(u). The imaginary part is non zero only in a finite range of frequency 
td-(T,n) < id < id + (T,n) that corresponds to the dispersion of the collective polarization 
erodes, regarded as a distribution of eigenvalues. Then dens,ty of state (DOS) is g.ven by 
D(uo) oc ida" (id) [25J. Fig. 2 illustrates our calculation at T = — 35°C. At low concentrations, 
the DOS is peaked at u (T), whereas the spectrum broadens progressively as the density n 
increases, indicating the drastic effect of the interactions. 

Since the squared eigenfrequencies of the collective modes are to be positive, the low edge 
Ld 2 _ (n, T) drives the stability of the system. Generalizing the Herzfeld criterion of polarization 
catastrophe, we define the critical density n c \ of the MIT as 

w_(T,n c i) = 0. (2) 

We have reported in Fig. 1 the calculated critical densities obtained with the use of Eq.(2). 
It indicates that the MIT occurs between 2 and 5 MPM depending on the temperature, 
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Fig. 4: Reduced isothermal compressibility of the electron gas as a function of concentration at 
T=— 70° C .The solid and the dashed curves correspond to the case of Na and Cs counterions, 
respectively. The arrow shows the locus of the dielectric catastrophe for the solvated electrons, 
indicating a miscibility gap between n c i and n C 2(Na) for Na, which does not occur in the case of 
Cs ions (see text). 

which is quite comparable to the experimental data. The evaluation of the original Herzfeld 
critical density, i.e. without taking into account the interactions, provides the MIT located 
at about 14 MPM at T = — 70°C. That shows how important the effect of the interactions is 
to correctly evaluate the MIT. Another consequence of the polarization catastrophe is that 
the low-frequency dielectric constant e(T,u>) diverges at n c i as it is experimentally observed 



in MAS 



27 



281 ] . Comparing the calculated data with the experimental ones, we find a good 
agreement between them (Fig. 3a). We also have calculated the real and the imaginary 
parts of the dielectric constant and evaluate the optical absorption coefficient A(u). Again 
the calculated concentration dependence of A(oj) at the locus u> max of its maximum, agrees 
well with the experimental data {29), [3^] at concentrations below n c \ (Fig. 3b). 

Coming at the issue from a different angle, let's consider the question : could the dipo- 
lar fluctuations lead to the phase separation experimentally observed in MAS ? Above the 
critical concentration n cl , the localized electrons are not stable. Hence, the behaviour of the 
system above n cl depends on thermodynamics of the metallic state. But the homogeneous 
electron gas is known to be unstable at sufficiently low densities due to occurrence of a 
negative compressibility. Therefore, if the polarization catastrophe occurs at a lower den- 
sity than this instability, it should provoke a phase separation. The MAS seems to be just 
this case. To reveal it, we have calculated the electronic part of the compressibility for the 



metallic state with the use of a modified model of stabilized jellium [3l| (see Methods) for 
several alkali metals. Fig. 4 shows the two curves obtained for Na and Cs respectively. In 
the case of Na, the compressibility Kp diverges at a critical density n C 2 ~ 6 MPM, whereas 
the dielectric catastrophe of the solvated electron state occurs at n c \ ~ 5 MPM. This is the 
origin of the miscibility gap and the associated phase separation : it exists a range of density 
n G [^ci,^c2] for which both states are unstable. Another consequence of this phenomenon 
is that a Mott mechanism for the MIT, which requires a stable electron gas at the critical 
density, appears impossible since the experimental MIT occurs at lower concentration than 
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n c2 . The second curve in Fig. 4 is for Cs. It is seen in that n c2 < n c i, contrary to the case of 
Na. No miscibility gap is thus expected in the case of Cs. This result is also coherent with 
the experimental facts. Although our estimations give upper and lower bounds for n c \ and 
n C 2 respectively, underestimating the instability range, they reveal a correct trend in the 
dependence on the size of ions, namely, a decrease of the instability range for heavier ions, 
due to scattering of delocalized electrons on ion cores. The later may decrease n c2 enough 
to destroy the miscibility gap. 

In conclusion, following the Herzfeld idea and using the hard-sphere models for quantum 
polarizable fluids, we have evaluated peculiarities of MIT in MAS, namely, the anomalies of 
dielectric response and concentration changes in the absorption maximum. Our estimations 
of the behaviour of the insulating and the metallic phases have revealed an instability range 
at low temperatures. Although the predicted miscibility gap is sufficiently smaller than 
the experimental one and is sensitive to variations of the model parameters (for example, 
10% decrease in eoo enhances twice this range), our calculations tell us that MAS system 
may deviate from the usual Mott scenario. The MIT seems dominated by the old Herzfeld 
mechanism. Our model is simple and does not take into account phenomena which may 
influence the transition such as disorder effects on the cavity formation of localized electrons, 
clusterization of electrons, influence of ionized states of metal atoms, and so on. These effects 
may be important in the vicinity of the transition, nevertheless simple estimations yield 
their energy scales sufficiently lower than that of the dipolar interactions between solvated 
electrons. Long-range nature of these interactions and low density of excess electrons in the 
solution support our game with the semi-analytical calculations. 

We also believe that the proposed scenario may be applied to some other systems. The 
reason of our optimism is the following. The behaviour of the quantum polarizable particles 
with respect to the homogeneous electron gas is controlled by two dimensionless parameters 
a a~ 3 and na 3 . The concentrations of electrons are to be low to provide the instability of 
the electron gas, and, hence, na 3 is to be small. But the parameter a o"~ 3 is also small 
in the case of ordinary polarizable fluids. For such systems, the present scenario remains 
unlikely. However in our case of self-trapped quantum particles, a^a^ 3 may be large enough 
to provoke a phase separation (it is about of 0.5 for MAS). Hence, our scenario is general 
and may take place in other systems, depending on the origin of the polarizable particles. 
For instance, in alkali metal-alkali halide solutions, solvated electrons, phase separation and 
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dielectric anomalies were experimentally observed [35]. Another example is the case of doped 
polar solids such as oxides where formation of large polarons occurs. For such materials 
aoa~ 3 ~ e^ 1 — e^ 1 and the origin of the unusual behaviour is essentially due to different 
scales of screening at various frequencies e s >> eoo. We have previously revealed it in the 
case of a Wigner crystal of polarons [7] which, from this point of view, behaves similarly 
like MAS. Finally, we would like to emphasize that our scenario opens a window in the 
general understanding of the MIT, facing a situation when the metallic and the insulating 
states are both unstable in a finite range of density. This stimulates the further challenging 
question. Can the dipolar interactions of localized quantum particles provoke a transition 



to a superconducting state ? Such a possibility was already discussed in [36[ , anc . 



perhaps, 



a careful analysis of former experiments on superconductivity in frozen MAS 37J, |38j could 
help to find the answer. 

Methods 

Evaluations of the effective polarizability 

In the case of interacting quantum Drude oscillators the problem is to evaluate the effec- 
tive polarizability a(u). Due to dipolar interactions between the oscillators, the polarizability 
a(u) is modified with respect to the non- interacting case and is given by the self-consistent 
equation : 

too/a (u) = eoo/ao (w) - 2E(a (u) /e^), (3) 
where the last term accounts the correlations between induced dipoles. Eq.(3) has been 



derived in 23] , we only modify it by taking into account the high-frequency screening by the 



solvent (the use of in Eq.(3)). Formally the quantity 3aE(a)/(3 equates with the dipolar 
part of the internal energy per particle of_a classical liquid of nonpolarizable particles with 
permanent dipole momentum (3a/ P) 1 ^ 2 
the function E(a) is known, all the physical properties of the system can be evaluate d by 
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26), where (3 is the inverse temperature. Once 



solution of Eq.(3). The simplest method to obtain E(a) is the Pade approximation [23|, |24j, 
which is an interpolation between the case of low and large polarizability a{uS). Adapting 
this method to our case, we get 

T7i / / \ / x I (na 3 )na(uj) 

E (a (u /eoo = ^^-7 — 3W r , 4 
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where Io(x) and h(x) are analytical functions depending on dimensionless density x = 
na 3 26j. Replacing Eq.(4) in Eq.(3) leads to a quadratic equation, which allows complete 
calculations of both the real and the imaginary part of a(u). Our input phenomenological 
parameters are 1) a — 2r c = 3.2 A [14] which remains fixed for all our calculations, 2) 
600=1. 76 33], and 3) uq(T) extracted from experimental data. This last parameter induces 
an implicit temperature dependence of a(u>). From [32], we take u} (T = — 70°C) = 0.9 
eV and du {T)/dT = -2.2 ■ 1(T 3 eV/K for higher temperature (-70°C < T < +70°C). 
The static dielectric constant in Fig. 3a is calculated by taking e s (T = — 70°C)=25, and 
de s /dT = -0.1 K- 1 [33]. 



Stability of delocalized states 

In a metallic state, the solution represents a plasma consisting of a degenerate electron 
gas strongly coupled with the ions dissolved in ammonia. Although a microscopic study 
of the system is still beyond possibilities of current methods, low temperatures and low 
metal concentrations of MAS simplify our analysis. We treat the influence of the solvent 
simply as a screening effect of the interactions, but we use different dielectric constants 
for the interacting electron gas and the ionic potential, because the latter are additionally 
screened by the orientational polarization of solvent molecules. We characterize the plasma 
by dimensionless parameters r s = (47m/3)~ 1/ ' 3 a^ 1 , T e = / 9e 2 /e 00 aBr s , and Tj = (3e 2 / 'e s aBr s , 
where as = h 2 /me 2 . We can express the change of the free energy caused by the dissolution of 
metal atoms in terms of the dimensionless parameters, and write the change Af per electron 
(or per metal atom) as the sum of the electron, the ion, and the electron-ion contributions 

A/(n) = A/ 6 (r 6 ) + Afi(Ti) + A/ ie (r e , r is r„ a,), (5) 

where is a parameter related with the short-range electron-ion pseudopotential, which 
takes into account deviations from Coulomb interactions between electrons and ions. Because 
T e ^> Tj, it may be checked that the electron gas gives the main contribution to A/(n), 
whereas the ionic contribution is only a correction. We can thus ignore thermal effects for 
simplified evaluations. We take the expression obtained for the stabilized jellium model as 



m 



3l|: 



A/(n) = - ^ + e c (n) + ?» + ^ ( 6 ) 
10 47reoo e.r. 
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where kp = {?>ix 2 n) l ^aB is the Fermi wave vector, and C m ~ —0.89774 is the Madelung 
constant in atomic unit (a.u.). The first, the second, and the third terms in Eq.(6) are the ki- 
netic, the exchange, and the correlation contributions to the total energy of the homogeneous 
electron gas with a positive jellium background. The last two terms are corrections which 
take into account the atomic nature of the cations. The only difference with the stabilized 



jellium model 311 ] . resides in the use of the dielectric constant in the relevant contributions. 



The local density approximation 34j was applied to calculate e c (n). The main difficulty is to 
evaluate a^. In simple metals, the parameter a« may be estimated as di m 2irR 2 /3, where R 
is the ion-core radius related with the atomic number of the ion. Advanced models treating 



smooth continuous pseudopotentials [3 11 ] yield numerical corrections to this trend. We apply 
(ii derived from [3l| and take the screening effect of solvent into account, i.e. a« = ai/e^. 
Importantly, we here use because a, represents the short-range part of the interaction 
between electrons and ions. As a result, we find at =10.8 and 26.1 a.u. for Na and Cs res- 
pectively. Numerical deviations of about 10 percent from these values do not change our 
results significantly. With the use of Eq. (6) we calculated the electronic part of the reduced 
isothermal compressibility Kp ~ [nd 2 (nj3Af)/dn 2 }~ 1 and found the concentration n c2 below 
which the electron gas is unstable (Fig. 4). 
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